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By applying the projection technique to the computation of excitation energies, 
we study the correlation effects on the band gap of conducting polymers. In the 
-vj . presence of an additional electron or hole, the correlation induces a polarization cloud 

around the additional particle, which forms a polaron. For the excitation energy of a 
►^ polaron, there is a competition between a loss of the correlation energy in the ground 

^SJ ' state and a gain of polarization energy. For the Hubbard interaction, the loss of 

^^ . correlation energy is dominant and correlations increase the band gap. However, for 

f^ , long-range interactions, the gain of polarization energy is dominant and correlations 

decrease the band gap. Screening the long-range interaction suppresses the gain of 
Qv I the polarization energy so that correlations again increase the band gap. A small 

dimerization is always favorable to the correlation effects. For frans-polyacetylene, 
we obtain the on-site repulsion U = 4.4 eV and the nearest-neighbor interaction 
V = 0.8 eV. The screening of vr electrons due to the polarizability of a electrons is 
quite strong. 
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I. INTRODUCTION 

In conducting polymers, such as trans-polyacetylene, the importance of electron-electron 
interactions has been widely accepted. Many physical phenomena, for example, the nonva- 
nishing negative spin density on alternate carbon atoms |jl|, the relative ordering between 
states 2^Ag and l^B^ 0, the optical absorption associated with neutral soliton 0, etc., 
should be interpreted by coping with electron-electron interactions, or more precisely, with 
electron correlations. But how strong are the correlations? This question has been argued 
for many years. A number of theories have been developed over the years to treat elec- 
tron correlation in polyacetylene, among them are the mean field and perturbation theories 
0J^J^, the valence bond analysis 0, the Monte Carlo calculation [|,3, Local Ansatz method 
|T^,[TT|, and the correlated-basis-function theory P^JTB[ ]. Comparison of these results with 



experiments has led to conflicting claims about the strength of electron correlations. 

The band gap is an important physical quantity which is strongly dependent on electron 
correlations. It can be determined by the optical absorption and for irans-polyacetylene, a 
value of 1.8 eV is found fl^- It is well known that calculated band gaps come out much too 
large when the independent-electron or self-consistent field (SCF) approximation is made. 

Recently, Konig and Stollhoff [|l^ performed an ab initio calculation for the ground state 
of polyacetylene. By fitting the dependence of the total energy on the dimerization, they 
were able to determine the on-site electron-electron interaction U = 11.5 eV as well as the 
effective interaction of electrons on nearest neighbor sites, V = 2A eV. When the energy 
gap is computed within SCF approximation a value of 6.9 eV is found, which is much larger 
than the observed one of 1.8 eV. The difference between these two values is apparently due 
to correlations. 



The projection technique |]I6|, developed recently by Fulde and his coworkers [|l^,|l3], has 
been successfully applied to both weakly and strongly correlated electron systems p!9| , |20 . 
In this paper, the technique is applied to the calculation of the excited states of conducting 
polymers. One advantage of the technique is that it works on r instead of k space. This 
enables one to generate with increasing accuracy the local correlation hole around an electron 
in the conduction band and a hole in the valence band and to interpret physically the 
different contributions to it. They are carefully studied both for an on-site Hubbard type 
interaction and for long-range interactions. The polarizability of a electrons enters the 
present calculations as an additional screening mechanism of the vr electrons. 

This paper is organized as follows. In the next section, we describe the model Hamiltonian 
which we use in this work. In Sec. Ill, the projection technique for excited states is described 
and the excitation energies are given. The competition between a loss of the correlation 
energy in the ground state and a gain of the polarization energy is shown. In Sec. IV, we 
give the correlation gap of conducting polymers under the scheme of the projection technique. 
In Sec.V, the numerical results are given, both for the Hubbard and long-range interactions. 
They are compared with the exact solution at the undimerized and the independent dimer 
limits. In the last section, we discuss the results we obtained and the screening mechanism 
of 71 electron interactions for trans-polyacetylene. In the appendixes, we give some details 
of evaluating and the gap corrections by second-order perturbation theory. 



II. MODEL HAMILTONIAN 

The vr-electron system of one-dimensional conducting polymers can be described by the 



Pariser-Parr-Pople (PPP) model pT 



H = HsSH + Hitter, (2.1) 

where, 

^SSH = -Z![^o - a{ui+i - ui)]{cl^i,cis + h.c.) + ^Kj^i^i+i - ui)^ (2.2) 

l,s I 

is the Su-Schrieffer-Heeger (SSH) Hamiltonian [Q, q^ and Qs are electron creation and 
annihilation operators at site / and with spin s{= a,/?), ui is the deviation of the /-th 
site from the equilibrium position (with equal distance between sites), the nearest-neighbor 
hopping integral between /- and / + 1-th sites has been taken as a linear function of the 
bond modulation ui+i — ui, and the last term in Eq.(2.2) is the elastic energy of this one- 
dimensional lattice. Furthermore, 

H,nter = ^U J^ PlPl + 1 E '^ll'PlPl' (2-3) 

I 1,1' 

is the electron-electron interaction, which includes an on-site Hubbard term U and off-site 
interactions Vw between electrons on site / and /', and pi = q^cia + CipCip — 1 the net charge 
density. The prime in the second summation term implies I ^ V . 

For the calculation of ground state, an effective short-range potential, i.e. Hubbard or 
extended Hubbard model (which includes nearest-neighbor repulsion) is commonly used to 
account for vr-electron interaction in polymers. However, when the system is away from the 
fuUy-shell electron distribution, long-range polarizations become important. So it is better 
to use a long-range potential for a correct description of vr-electron interactions in excited 
states. The interaction range depends on the screening by the polarizability of a electrons 
in the system. Various empirical formulas |^ for this long-range potential of vr electrons 



have been worked out. Here we use Ohno's formula 12^J2T| for the interaction VJ//, which is 
given by 

Viv = , ^ (2.4) 

v/l + {ru'/roy 

Here, tq ^ 1.29A, r;// is the distance between sites / and /' in unit of A and U = 11.13 eV. 
Roughly speaking, for polyacetylene, the carbon bond angle is 120° and average distance 
between the nearest-neighbor sites is 1.4A. 

For a dimerized lattice {ui = ±m), the SSH Hamiltonian can be diagonalized exactly. 
Within Hartree-Fock approximation, we can obtain a self-consistent solution by including 
the exchange parts of the interaction (2.3). However, instead of a self-consistent HF Hamil- 
tonian, we simply start from the following effective single-particle Hamiltonian, 

Ho = -tE[l ± ^](4+i.Q. + h.c), (2.5) 



where, 4t is the energy band-width and z is the gap order parameter, + for long bonds 
and — for short bonds. Physically, there are two factors to influence the values of these 
two effective quantities t and z, besides the bare single-particle contribution, which gives to 
and 2au/tQ respectively. One is the exchange-interaction contribution, which results from 
off-site interactions. For the extended Hubbard model, this exchange-interaction gives 

t = to + Vm, (2.6a) 

2au + V6m 
z = , 2.6b 

where rn and Sm are defined by (c;\,q+is) = fn±5m [^. For a long-range interaction, we only 
consider an effective change on the nearest-neighbor hopping instead of explicitly including 
these hopping terms beyond the nearest-neighbors |Ty]. Other is the correlation contribution 



ri| , which reduces the band-width and increases the gap order parameter. While the treat- 
ment of exchange contributions is Hartree-Fock approximation, the inclusion of correlation 
contributions implies that we have gone beyond the mean field theory. Mathematically, t 
and z are determined by minimization of the total (both mean-field and correlation) en- 
ergy expectation value of the Hamiltonian (2.1) in the ground state. It has been known 

111 that the correlation contribution to the parameter z becomes important and makes the 
total energy much lower in the interaction region {U > 2t), and with this contribution the 
calculation using local ansatz is valid until intermediate interaction {U ~ 4t). 

In the Hamiltonian (2.5) we divide the operator cis into two parts, i.e., the electron 
operator in conduction band ais and the hole operator in valence band bis, 



where, 



cis = ais + h\s, (2.7) 



ais = ^=^e*'''[2sin0fc + {-lY cos^kWs, (2.8a) 

4 = ^ E e^''[cos 0. + ^(-1)' sin 0,]6L,„ (2.8b) 



s = —s, the summation is over a reduced Brillouin zone (— 7r/2, 7r/2], and the lattice constant 
has been taken as unit. The Hamiltonian (2.5) becomes 



Ho = 2tJ2<k){alaks + blhs - 1), (2.9) 

k,s 

and the electron energy spectrum (in unit of 2t) 



e{k) = Vcos2 k + z^ sin^ k. (2.10) 

The transformation angle 0^ is defined by 



For a half-filling case, the ground state |0) is the vacuum of the electron operator aks and 
hole operator b^g, so that it is also the vacuum of operators a/s and big. 
We define the correlation functions Pw and Qw |]T0[ as follows 



Pw = {clcvs) = ik-sbU), (2.12a) 

Qw = {ciscl,) = {aisal), (2.12b) 

and the anticommutators for the operators a and b are 

{ais,al^,} = 6ss'Qii', (2.13a) 

{bi,,bl^,}=6ss'Pw. (2.13b) 

Explicit expressions for these functions are given in Appendix A. We will need them in the 
calculations of the following sections. 

Substituting Eq.(2.7) into the interaction (2.3), we have the two-particle part as follow 



where, 



and 



H, = J2VmOm, (2.14) 



O^ = 0(^) + (0^2) + h.c) + (Oi^) + h.c), (2.15) 



OU = T,lss'{{0'ls(^l+ms'(^l+ms'ais + blsbl+ms'bl+ms'hs) 

-{O'lsh+ms'bl+ms'O'ls + ai^^gbi^,bis'ai+rns) 

+ («L4^«+ms'az+ms' + aUmsbUmsbis'ais')} (2.16a) 

is the interaction between electron-electron, hole-hole and electron-hole; 

OL'^ = E ^Ulms'bl^A (2. 16b) 

Iss' 

is the spontaneous creation of two electron-hole pairs and 

^m — Zeiss' \(^ls^l+ms'"l+ms'^'-^ ~^ '^l+ms'^ls'^isi'^l+ms 

+b\sh\+rr,s'0^l+ms'^ls + b\+^sb\s' a\s,h+ms) (2.16c) 

is the creation of an electron-hole pair through the scattering of an electron or a hole. We 
have taken the translation invariance for the interaction coefficients, i.e., Vm = VJ^^+m for 
m 7^ and Vq = U /2. The one-particle part of interaction (2.3) has been absorbed into the 
effective Hamiltonian (2.5) as the exchange contribution. 



III. PROJECTION TECHNIQUE 

We consider the case an electron or hole is added to the half-filled system. The excitation 
energy is defined by 

e{k) = E^{k) - Eq (3.1) 

and is contained in the one-particle correlation function 

i?.(fc,r) = (^|c..e--(^-^°)cL|^). (3.2) 

Here, \g) is the ground state of the half-filled system, Eq is the corresponding energy, c|.^ 
creates an electron {a\g) or a hole {h\g), and Ei{k) is the energy for the system with the 
additional particle. Within the quasiparticle approximation, the correlation function reduces 
to 

The quasiparticle energy e{k) shows up as a pole of its Laplace transformation Rs{k,x), 
which has been shown [0] to be 

R.ik,x) = {cln\——^——cln). (3.3) 

X — [Ijo -\- til) 

In the above, we have used a bilinear form in Liouville space 

{A\B) = {0\A^B\Or, (3.4) 



the index c expresses the cumulant production |T^, which is defined in Appendix B. The 
Liouville operator Lq acting on any operator A gives 

LoA=[Ho,A]_, (3.5) 

and 

n = lim(l + J H,), (3.6) 

^•^0 X - {Hi + Lq) 



which transforms the ground state of Hq (|0)) into the exact ground state {\g)) 

Now, we divide the Liouville space R into a relevant part Rq and an irrelevant part -Ri. 
Let Rq be spanned by a set {A,y} of elements \A,^). Then the projector |]16| 



P = j:\A,)x-JiA,\ (3.7) 






projects onto the space Rq, and Xfiu = {A^\A^). 
Within this relevant space Rq, we have 



Rs{k,x) = J2{cln\A,)x-jR.{k,x), (3.8) 



fl,!/ 



where R^ik^x) denotes the set of functions 

Mk,x) = (A^l ^-^^Pcln). (3.9) 

Then we have the equation 

Y^ixxuv - ^•yv)x^,^Rt^{k, x) = {A^\Pci,Q) (3.10) 

with 

uj^, = {A^\{Lo + H,)A,). (3.11) 

The excitation energies e{k) are given by the poles of Rs{k, x) or Ru^k, x), so are the solutions 
of the following equation 

det{xXf^u - uj^i,} = 0, (3.12) 

if we choose the elements A^, which span the relevant space Rq, so that (A^lcl^Q,) ^ 0. 

In the half-filled ground state, the electron distribution is fully-shell. The correlation 
interaction (2.14) would give rise to spontaneous excitations of electron- hole pairs. These 
spontaneous excitations lower the total energy. In the presence of an additional electron or 
hole, the correlation causes two processes: one is the blocking of the spontaneous excitations 
of electron-hole pairs if these excitations are associated with the extra particle; other is the 
excitation of electron-hole pairs through the scattering of the extra particle. The former 
(which we use operator S"^ to describe) gives rise to a loss of the correlation energy in 
ground state so that the excitation energy is increased. And the latter (which we use 
S'l^ to describe) is a polarization process, which decreases the excitation energy by the 
gain of the polarization energy. These two processes induce a polarization cloud around 
the extra particle so that an electron polaron or a hole polaron is formed. Here we set 
{A^} = {c|.^} © {_S'^ii(^^!^ © {cL'S'^}) where {5*^} are the operators, which do not commute 
with the extra particle creation operator c|,^ while {S''^} commute with the operator c^^. 
Then we can rewrite the excitation energy as two parts 

£(A;) = 2t-e(A;)+£"°"(A;), (3.13) 

i.e., the mean-field energy spectrum and correlation contribution. The latter is combined 
by two parts, too, 

e^°"(A;) =£'^(/t)+e''(A;), (3.14) 

where 

e-{k) = -E^m(^)(4|^1'5;4) (3.15a) 

is the gain of the energy from the polarization process with the extra particle, and 

e\k) = Y.^,{k){c{,\H,c{,Sl) (3.15b) 



is the loss of the correlation energy in the ground state due to the extra particle. 

When the mutual influence of the two types correlation is neglected by set 
i^Ji'^ksl-^i'^ks^u) — 0) which we believe is small and is negligible for a non-strong interaction 
(see Appendix C), then we can write down 

n,ik) = E^;.'(^)(^:4l^i4), (3.16a) 

V,{k) =T.C;.'iv)icls:j\H,cl). (3.16b) 

V 

The coefficients are the inversions of matrices C(7r) and C{ri), which are given by 

C,M = iS;cl\[Lo + H,- Eik)]S:cl), (3.17a) 

C,M = -{clsmo + H,- eikM,S:j). (3.17b) 



When the mutual influences of the correlated excitations and the polarized scatterings are 
taken into account, the expressions for vr^(fc) and i]^{k) become somewhat more complex. 

From the expressions (3.17), we notice that both the loss of the correlation energies and 
the gain of polarization energies are dependent |]T6[ of the exact excitation energy e{k). This 



dependence comes from the non-zero of susceptibilities Xij.u- It makes the excited energies 



must be obtained self-consistently and results in a narrowing of energy bands |2^ . The more 
important is that, we will see in Sec.V, the dependence ensures that its results are correct 
for intermediate interaction strength. 

IV. CORRELATION GAP 

A band gap in a semiconductor is the energy it costs to move an electron on the top 
of valence band to the bottom of conduction band. In an electron-hole picture it is the 
minimum energy to create an electron-hole pair. Here, we deflne the correlation gap Eqg as 
a contribution of the correlation effect, i.e., it is the difference between the exact band gap 
and the mean-fleld band gap 

Egg = -Egap - -Egap. (4.1) 

For our case where there exists the symmetry of electrons and holes, the excitation energies 
of an electron polaron and a hole polaron are same. So the correlation gap should be twice 
of the excitation energies of the polaron ^^''''(fco) {ko = vr/2), if the interaction between 
the electron polaron and hole polaron — the excitonic effect — is neglected. We leave this 
excitonic effect in this system as a separate study. 

For the calculation of this excitation energy, we should flrst choose the polarized scat- 
tering operator S^l and the correlated excitation operator S^. When the interaction is not 
strong, we can choose 

S:, = '^T.4.blal,aj,^ (4.2a) 

a, a' 

and 



^1^ = IE«144.'^1-'- (4-2b) 

These operators span a relevant space Rq for the case an electron is added into the conduction 
band. The physical meaning of the two kinds of operators has been discussed in the previous 
section. 

The choice of the above operators means that we do not break the spin symmetry. It is 
correct for a dimerized system, where mean-field ground state is of bond order waves (BOW) 
or charge density wave (CDW), without spin density waves (SDW). However, for a uniform 
lattice — the undimerized limit — the SDW would play a non-negligible role. It leads the 
calculation using the above operators underestimate the energy gap. Furthermore, we have 
only considered two electron-hole pairs excitations, which is enough only for a non-strong 



interaction. As shown in the calculation for ground state, it works until f/ ~ 4t ||Tl|. This 
is enough for the polymers like polyacetylene, where U ~ 4to- We will see it also holds 
for excited states in the following section. However, for U > At, the system will undergo a 
spin-Peierls transition, there the spin-spin interaction is dominant. This is not the regime 
we interest in this paper. 

With the translation invariance, we reduce the dimension of relevance space Rq by using 
the following operators, 

S: = T.Sh^ci, (4.3a) 

i 
i 

Now we write down the following four kinds of matrices: 1. susceptibility- matrix 
(S^aislSd'aL), ('^>Ll'^>L); 2. projection matrix {aljHiS^al^) , (4ji:^i^>L); 3. hop- 
ping matrix {S^aljHoS2,al,), (5'>L|ifo5'>L) and 4. interaction matrix (S'Ja^lifiS'JaL), 
(SjalglHiSj/al^) . They should be evolved before the calculation of excitation energies. 
These matrices are given in Appendix B. 

V. NUMERICAL RESULTS 

A. Hubbard interaction 

We consider Hubbard interaction in this subsection based on the general formula estab- 
lished in the previous two sections. For such a short-range interaction, we first choose only 
on-site local operators SJ and 5*^ with d = as the elements of set {^/i}- With this choice, 
we have the expressions for excitation energies 

where, the coefficients a and b are as follows 



a^ = Y.RiPj', (5.2a) 

I 

b^ = 3j2EiRiP^/a^; (5.2b) 

I 

and 

a, = Y.i-lYRiPJ', (5.3a) 

br, = Ate{ko) + 3 J2{-iyEiRiP^/ar,. (5.3b) 

I 

The functions Pi, Ei and Ri are defined in Appendix A. The coefficients 6,^ and 6^ are 
independent of Hubbard strength U because the interaction matrix vanishes here. 

The equation (5.1) combined with the equations (3.14) and (3.15) determines the ex- 
citation energy e^°"{kQ) as well as the correlation gap Eqg self-consistently. For trans- 



poly acetylene, the band gap has been measured as 1.8 eV [0]. For Hubbard interaction, 
we will see the correlation gap Ecg is always positive, i.e., the band gap is enhanced by the 
Hubbard correlation. At ;z = 0.15, which opens a band gap of 1.5 eV in the single-particle en- 
ergy spectrum if t = 2.5 eV as taken by SSH |^, we have a-,, = 0.055, a^ = 0.196, 6^ = 4.34t, 
and brj = 4.43t. The correlation gap Ecg obtained is plotted as a function of U in Fig.l as a 
dashed line, while the result obtained with off-site operators d up to 10 is shown as a solid 
line. As to make a comparison, we have shown the second-order perturbation result as a 
dotted line and the Monte Carlo results P] as points with error bars in Fig.l. 

First let's see the limit of small U, where the correlation gap is of the form nU'^/t. The 
result with on-site operators gives n = 0.0677, which is about 90% of the value n = 0.0755 
calculated from the second-order perturbation theory (see Appendix C). Here we should 
mention that the second-order perturbation theory for the correlation gap is exact in the 
limit of small U for the case a finite dimerization has opened a gap. It is different with 



the uniform lattice {z = 0) case where the band gap is purely Coulomb gap [26| and the 
perturbation theory doesn't work. Our result can be improved greatly by include of off-site 
operators, for d up to 10, we obtain k = 0.0723, which reaches 96%. From Fig.l, we can 
see that the off-site operators become more important when the interaction strength U goes 
larger, which is easy to be understood. 

At z = 0, the undimerized limit, the exact gap is available from the paper of Lieb and 



Wu 26 
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EcG = U-4t + 8 V(-l)"[inf/ - it' + WUY% (5-4) 



n=l 



which is shown in Fig.l as a dot-dashed line. Compared with this exact solution, the 
correlation gap a.t z = 0.15 obtained by our projection technique is always larger even only 
with the on-site operators. It means the correlation gap is enhanced by a finite dimerization. 
However, the perturbation theory shows us that there is a critical interaction strength Uc, 
when U is larger than Uc, the dimerization reduces the correlation gap while it enhances the 
correlation gap when U is smaller than Uc- Obviously this is an artifact. The fact tells us it 
is important that the dependence of both the correlation energy and the polarization energy 



10 



on the exact excited energy |2^ because it is the dependence which raises the correlation 
gap. 

In Fig. 2, we show the dependence of the correlation gap on the dimerization z for different 
interaction strengths [/. At very small dimerization 2;, as declared previously, the result we 
obtained here underestimates the correlation gap because the choice of the operator set 
(4.3), which doesn't break the spin symmetry. Comparing with the exact solution at 2; = 0, 
we can see that this region is very small for small f/, and it is about 0.1 for \J = 4t. 
The enhancement of the correlation gap by a small dimerization is shown clearly and the 
correlation gap reaches its maximum around z = 0.3. For a larger z, the correlation gap 
decreases monotonously. This behavior is easy to be understood by compared with the 



correlation energy in the ground state [|T0[, which is decreased by a small dimerization and 
after reach its minimum value it goes up monotonously with the increase of the dimerization. 
The reason is that for Hubbard interaction, the loss of correlation energy is dominant and 
the polarization energy is small. 

At 2; = 1, the independent dimer limit, the exact correlation gap can be obtained easily, 

EcG = {U' + 64t'Y/'-8t, (5.5) 

which has been indicated in Fig. 2 as well as the exact solution (5.4) at 2; = 0. At the 
independent dimer limit, the long-range operators are not important while the high-order 
excitations become important for a strong interaction. Since the calculation only with low- 
order excitations gives a value larger than the exact result, the high-order excitations should 
reduce the correlation gap, which is contrary to the effect of long-range operators. So we 
expect that there are some cancellations between the longer-range excitations and the high- 
order excitations. 



In the end of this subsection, we mention the work of Sun et al. [12], which used a screened 



Coulomb interaction and found the optical gap is reduced by the electron interaction if the 
screening is very strong. This result is not contrary to ours because as stated in their paper 
this reduction is due to the non-diagonal bond-charge repulsion, which we don't include in 
this paper. 

B. Long-range interaction 

For the long-range interaction, we take the Ohno formula (2.4) to parameterize the 
interaction coefficient Vw, which is dependent of the distance between the site / and V . 
With a dimerization u, the change of bond lengths is 4m/v^ and the bond alternation is 
formed. Within the mean-filed approximation and only nearest-neighbor-site interaction, 
the relation of the dimerization u and the gap parameter z is shown in Eq.(2.6). With the 
longer-range interaction as well as the effect of correlation, the dependence of u and z should 
be determined numerically as discussed in Sec. II. 

Due to the existence of the long-range interaction, the relevant space Rq only with the 
on-site operators {d = 0) is too small to give a qualitative correct result, the inclusion 
of long-range operators is necessary because the polarization process is of a much longer 
correlation length than that of correlated excitations. With the relevant space Rq spanned 
by the operators (4.3) with d up to 10, the correlation gap is shown in Fig. 3 (a). We can 



11 



see that it is the long-range interaction which makes the gain of the polarization energy 
becomes dominant so that the correlation gap is negative. This result means the long-range 
correlation reduces the mean-field energy gap, which is unphysically large from the ab initio 



self-consistent filed (SCF) calculation [|T5l. This reduction becomes larger with the increase 
of a small dimerization. At about z = 0.1, the correlation gap has reached its minimum and 
then increases with the increase of the dimerization. 

At z = 0.2, the reduction of the band gap is about 3.4 eV. For trans-polyacetylene, the 
SCF gap is 6.9 eV |T^, among which the SCF dimerization contributes 1.05 eV and the 



exchange interaction contributes the rest. So the Coulomb gap arose from the electron- 
electron interaction is 2.45 eV. The deviation comes from the model interaction (2.3) which 
doesn't include the polarizability of a electrons in the system. The presence of a electrons 
would be favorable to the vr electron excitations so that it should reduce the band gap 
further. 

In a simple way, the polarizability of a electrons could be considered as a screening of 



the TT-electron interactions, i.e., the Ohno interaction is screened as follows 10 



V,,' 
Vu> = — . (5.6) 

Co 

At this moment, the on-site Hubbard interaction U is kept unchanged. With the screen for 
eo > 2, the correlation gap becomes positive again. It means that the screen suppresses the 
polarization process so that the gain of polarization energy is reduced and the loss of the 
correlation energy is again dominant. This is same as the correlation gap for a short-range 
(Hubbard) interaction qualitatively. Fig.3(b) shows the correlation gap for the screened 
interaction with eo = 3. 

The screening reduces the tt electron interaction range so that the correlation increases 
the gap. But at the same time, the screening reduces the exchange contribution, too, by a 
factor of I/cq. In fact the screen decreases the energy gap, i.e., the contribution of the a 
electrons tends to reduce the band gap. 



VI. DISCUSSIONS AND SUMMARY 



First let's see the ab initio work of Konig and StoUhoff [0. At SCF level, it gives the 
SCF energy gap 6.9 eV at the equilibrium dimerization u = 3.28 pm. The equilibrium 
dimerization u decreases to 2.52 pm when electron correlations are included. By fitting 
the ground state energy, they obtained the model parameters a = 40 meV/pm, K = 3.9 
meV/pm^ for to = 2.5 eV. With the above parameters, we obtain the SCF gap of 7.0 eV 
for the long-range interaction (2.4). Among it, the dimerization contributes 8au = 1.05 eV 
while the exchange contribution 2J2i{—^)^^y2i+iP2i+i is 5.95 eV. The agreement with ab 
initio calculation shows that the PPP model (2.1) describes the vr-electrons in conducting 
polymers quite well. 

With electron correlations, the bond alternation is reduced [^. Then the dimerization 



contribution 8au to the band gap becomes 0.8 eV and the exchange interaction contribution 
to the band gap decreases 5.83 eV. With the screen (5.6), the exchange contribution to the 
gap is 2.49 eV for eo = 2 and is 1.52 eV for eo = 3. For the screen with a fixed on-site 
repulsion U (= 11.13 eV), the correlation gap is 4.58 eV ai e^ = oo {i.e., only Hubbard 
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on-site interaction and no exchange contribution), and —3.39 eV at eo = 1 {i.e., no screen). 
That gives the band gap 2.9 eV at the hmit of unscreening interactions (eo = 1) and 5.4 
eV at the hmit of completed screening interactions (eo = oo). Obviously in this way, the 
contribution of a electrons to the band gap is contrary to the physical fact we expect, i.e., 
the polarizability of a electrons should reduce the band gap. 

We need to consider a more realistic situation of conducting polymers: the on-site vr 
electron interaction is not completely local, i.e., it has a nonlocal part. So the on-site 
repulsion U would be screened by the polarizability of a electrons from its "bare" value 
(11.13 eV). We take the form U = Uq + V, where Uq (= 3.6 eV) is the "net" on-site 
repulsion and V is the nearest-neighbor interaction, which is screened by a electrons as in 
Eq.(5.6). In this way, at the limit of eo = C)0, the interaction is the Hubbard model with Uo, 
the corresponding correlation gap is 0.42 eV, and then the band gap is 1.22 eV. That shows 
a reasonable reduction from the polarizability of a electrons. To get the band gap 1.8 eV, we 
have to set eo = 9, a quite strong screening to vr-electrons, then the exchange contribution 
is 0.43 eV and the correlation gap is 0.57 eV. With the screen, the on-site repulsion U = 4.4 
eV and the nearest-neighbor interaction V = 0.8 eV, the interaction range is very short. 

The electron-electron interaction parameters U and V we determined from the energy 
gap are quite different with the values obtained by fitting the groundstate energy curve of 
the ab initio calculation on dimerization. In fact it has appeared for many years that the 
discrepancy of the strength of electron interactions among the theoretical results based on 
different physical quantities. Some of the results are tabulated in Table I. It seems that it is 
not appropriate to talk about the strength of electron interactions in polymers based on such 
a simple model. The physical reason behind it is that the effective vr electron interaction is 
strongly affected by the polarizability of a electrons, which play a different role for different 
physical processes. The interaction depends on physical processes. From Table I, we can 
see that there are mainly two kinds of physical processes for the interaction, one is in the 
ground state, other is in excited states. On the ground state, the interaction is about 10 eV 
]n|,[15| while it is about 4 eV [|,|H9],|0HT1 on the excited states. 



Finally, we summarize this work as follows. By using the projection technique for excited 
states, we have studied the correlation effects on the band gap of conducting polymers. In the 
presence of an extra particle (electron or hole), the correlation induces a polarization cloud 
around the extra particle, which forms a polaron. For the excitation energy of a polaron, 
there is a competition between a loss of the correlation energy in the ground state and a 
gain of the polarization energy. For a long-range interaction, the gain of polarization energy 
is dominant and the correlation decreases the band gap. A screening by the polarizability of 
a electrons reduces the interaction range and suppresses the gain of the polarization energy 
due to the correlation. The loss of correlation energy becomes dominant and the correlation 
increases the band gap for the short-range or Hubbard interaction. A small dimerization 
leads to a further enhancement (reduction) of the correlation gap for the Hubbard (long- 
range) interaction. For irans-polyacetylene, we obtain the on-site repulsion U = 4.4 eV and 
the nearest-neighbor interaction V = 0.8 eV by fitting the band gap to 1.8 eV. The screening 
of vr electrons due to the polarizability of a electrons is quite strong. 
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APPENDIX A: CORRELATION FUNCTIONS 

The correlation functions Pw and Qw defined in Eq.(2.12) have the relation 

Pw + Qw = Si,. (Al) 

For even value of m and n, 



") 



here 



J^m,n -* m+l,n+l 2 "^i^' (^A/aj 



1 /■'^/2 cosA;cos(2/ + 1)A; + 2;sinA;sin(2/ + 1)A; 

P21+1 = - / dk / 9, , 9 • 2 , • ^^^) 

TT Jo Vcos^ k + z^ sm k 



Vcos^ k + z^ sin^ k 
The energy correlation function is defined as 



= {brsHobU)-{H,){brsbU) (A4) 



and then 



El, = (cllHocl) = (allHoal) 

= i-iy^'Eu,, (A5) 

where, we have use the cumulant production. Similarly as did for Pur, for even values of m 
and n we have 

^mjU -^m+l,n+l ^m—n) yi\\yaj 

Em+l,n = En,m+1 = -C^m-n+1; (AoDJ 



here. 



and 



E2I = |t[(l + Z){P21+1 + P^2l+l) + (1 - Z){P21^1 + P^2l-l)], (A7) 

E21+1 = it[(l + z)5i,o + (1 - z)5i,^i]. (A8) 
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The polarized correlation function Rii'{k) is defined by 

{al\alaisal) = Ru'{k)/N, (A9) 

N is the number of sites. Similarly, we have 

ibllblksbl) = i-ir' Ru'{k)/N, (AlO) 

Ru'ik) = R:,,{k). (All) 

For even numbers of m and n, 

Rm,n\k) = Rm+l,n+l{kj = Rm-n{k), (Alia) 

Rm+l,n{k) = Rm^n+1 (k) , (Al2b) 

where, 

R2i{k) = e^'"', (Al3a) 

i?2m('t) = -e'''^''^'^"'''^'- (Al3b) 

At ko = 7r/2, R21 = —R21+1 = (—1)', independent of the dimerization z. 

APPENDIX B: COMPUTATIONS OF THE FOUR KINDS OF MATRICES 



The definition |16| of cumulant productions is as follow 



(0|A? ...A:\Or=^...^ ln(0| n e^'^lO)U,=...=,„=o. (Bl) 

Then we have the following properties: 

{0\AB\Oy = {0\AB\0) - (0|v4|0)(0|5|0), (B2a) 

{0\ABC\oy = {0\ABC\0) - {0\A\0){0\BC\oy - {0\B\0){0\AC\oy 

-{0\C\0){0\AB\Oy - {0\A\0){0\B\0){0\C\0), (B2b) 



In terms of the properties of cumulants and the correlation functions in Appendix A, we 
compute the four kinds of matrices as follows. 

1. Susceptibility matrix: 

iSIjakslSi'j'O-ks) = -jN^ii'^n'C^Qii'Qjj' ~ Qij'Qji'), (B3a) 

(SijakslSi'j'Ctks) = -4^Qii'Rjj'{2Pii'Pjj' - Pij'Pji') 

+ (,-^-) + (,'-/) + (,,,'-^-,/). (B3b) 

2. Projection matrix: 

(allOmSJ^al) = 2j2{SZ^^al\S:^al) + (m ^ -m), (B4a) 

I 

{^ksl^rni^ijO-ks) = '^Z^Wll+m'^ksl^ij'^ks}- (B4b) 

I 
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3. Hopping matrix: 

{Sij'^ks\^0'^i'j''^ks) = —j^{Rii'Pjj'{'^Eii/Qjj> + 2Qii>Ejjr — Qij'Eji' — Eij>Qji>) 

+Rii'Ejj/{2Qii'Qjj' — Qij'Qji')}, (B5a) 

Wij'^ks\Ho^i'j''^ks> — '^^^k{jij<^ks\^i'j''^ks) 

+RiiiEjji{APiiiPjji — 2Pij'Pjii)] 
+ {t ^ j) + it' ^ f) + {t, i' ^ 3,f)}. (B5b) 

4. Interaction matrix: 

\Sij^ks\OmSiijiaj^g) = —-^Rii' J2l {Pjj'Qli'Ql+mj' {'^QliQl+mj — QljQl+mi) 

'-^jj'^lj'^l+mi'y^^lj^l+mi '^li'^l+mj ) 
—'^Pl+mjPl+mj'Qli'i'^QliQjj' — QljQij') 
~^Pl+mjPl+mj'Qlj'\^QljQii' ~ QliQji') 
+4:PljPl+mj'Ql+mj'{'^QljQii' — QliQji') 
~^PljPl+rnj'Ql+mi'\2QljQij' ^ QliQjj') 

+ {m ^ —Tn)}, (B6a) 

i^ij'^ks\^mS,i,j,ai^^) = J^Y.1 {{RliQw + QliRli')Ql+mjQl+mj'{'^Pii'Pjj' — Pij'Pji') 

-\-Qii' Rjj' Pw Pl+rnj'i/'PliPl+mj ~ PljPl+mi) 

[\'^l+mi'^l+mi' ' ^l+miJ^l+mi' )^ jj' ~r ^l+mi^l+mi' -^jj'\ ^ 
x[Pli'{'^PliPjj' — PljPij') + Plj'i'^PljPii' — PliPji')] 
+ [{RliQl+mi' + QliRl+mi')Qjj' + QliQl+mi'Rjj'] X 

x[2-P«+mi'(2-PiJ-Pjj' — PljPij') '~ Pl+mj' {'^PliPji' — PljPii')] 

+ {m ^ —m)} 
+ (, - j) + (,'-/) + (,,,' -_^;/). (B6b) 

The matrices with the translational invariance are easily got from the above matrices. 

APPENDIX C: SECOND-ORDER PERTURBATION THEORY 

By use of the second-order perturbation theory, the correlation energy for the state with 
an extra electron in the half-filled ground state is 

m ^m ^k 

With the equation (2.14), we have the Hubbard interaction as 

H^ = W{Oi^^ + of + Oj^^^ + oS') + OS')^}, (C2) 
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and 

^i«L|0) = |f/{Of + O(=^)}aL|0), (C3) 

where, the first term gives the creation of two electron-hole pairs and other term gives the 
creation of an electron-hole pair through the scattering of the extra electron. So that there 
are two types of intermediate states, which have contributions to the correlation energy 
(CI). The contribution from the intermediate state 

1"^) = alis^LAsCsaLlO) (C4a) 

gives the correlation energy of the half-filled ground state and the loss of correlation energy 
due to the extra electron. The contribution from the intermediate state 

gives the gain of the polarization energy due to the polarized process after the additional 
electron was introduced. Then we rewrite the correlation energy in Eq.(Cl) as 

E^ik) = E-- + E.ik) + E^{k), (C5) 

where, E"^"^"^ is the correlation energy of the ground state, the loss of correlation energy in 
the ground state 

f/2 1 

^^(^o) = TTTT^ y" 7 ■ ■ ■ ^ (eoeie2e3 + C0C1C2C3 + Z0Z1Z2Z3 

4^^ yfctfca (^0 + ei + £2 + e3)eoeie2e3 

-CoZiZ2C^ - Z0C1C2Z3 + 2C0C1Z2Z3 + 2zoZiC2Cs), (C6a) 

and the gain of the polarization energy 

f/2 1 
En{ko) = --TTTH Yl 7 \ : ^ (eoeie2e3 - C0C1C2C3 - Z0Z1Z2Z3 

4^^ fcTfca (^1 + ^2 + £3 - 60)60616263 

+C0Z1Z2C3 + Z0C1C2Z3 - 2C0C1Z2Z3 - 2Z0Z1C2C3). (C6b) 

Here, we have defined k^ = ki + k2 — ko, ci = cos/cj, Zi = zsinki and 6^ = {cf + zfY^"^. So 
the correlation gap is 

EcG = 2Er,i7T/2) + 2E^{7T/2) 

= K{z)U^/t, (C7) 



where 



/ X _ 1 \-^ (- 616263 + Z1Z2Z3 — CiC2^3 + 2Z1C2C3 
'^y^) - ^^ 2^^ 1 (61 + 62 + 63 + ^)6i6263 

_ 616263 — ^i^2^3 + CiC22;3 — 2^:1 C2C3 

(ei + 62 + 63 - 2;)6ie263 



} (CJ 



and ks = ki + k2 — 'n 12. The coefficient k can be evaluated numerically. 

In the end of this Appendix, we mention that, in the second-order perturbation theory, 
there is no mutual influence between the correlated excitations and the polarized scatterings, 
which we assume is small and is negligible until intermediate interaction strength in Sec. III. 
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FIGURES 

FIG. 1. The correlation gap Ecg for the Hubbard interaction. The sohd hne is the result 
obtained with long-range operators (4.3) {d up to 10), the dashed line is the result with on-site 
operators (4.3) {d = 0), and the dotted line is the perturbation result at z = 0.15. And the 
dot-dashed line is the exact correlation gap at the undimerized limit z = (Ref. [^). The points 
with error bars are the Monte Carlo result (Ref. (P). 

FIG. 2. The correlation gap Eqg for the Hubbard interaction. The dots are got by the projec- 
tion technique. The lines guide the eye. The circles are exact solutions. 

FIG. 3. The correlation gap Eqg for a long-range interaction: (a) without screen eo = 1; (b) 
with screen eg = 3. 
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TABLES 



TABLE L The strength of electron interactions obtained by fitting to different physical quan- 
tities. 



U (eV) 


V (eV) 


Physical quantities 


11.5 


2.4 


The groundstate energy curve E{u)'^ 


7-9 




The dimerization u and the screened force constant*^ 


3-4 




The spin density of a neutral soliton^ 


3.7 


0.4 


The resonant Raman spectra*^ 


4.4 


0.8 


The band gap Eg'^ 



3.7-6.0 



0.0-2.5 



The band width W , u, Eq and the LO mode z^lo for N = 8, 10^ 



^Ref. H. 
^Ref. 111. 



^Ref. II . 
<^Ref. |29|. 
^Present work. 
^Ref. 
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